Kinetic roughening- like transition with finite nucleation barrier 
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Abstract 

Recent observations of the growth of protein crystals have identified two different growth regimes. 
At low supersaturation, the surface of the crystal is smooth and increasing in size due to the 
nucleation of steps at defects and the subsequent growth of the steps. At high supersaturation, 
nucleation occurs at many places simultaneously, the crystal surface becomes rough and the growth 
velocity increases more rapidly with increasing supersaturation than in the smooth regime. Kinetic 
roughening transitions are typically assumed to be due to the vanishing of the barrier for two- 
dimension nucleation on the surface of the crystal. We show here, by means of both analytic mean 
field models and kinetic Monte Carlo simulations that a transition between different growth modes 
reminiscent of kinetic roughening can also arise as a kinetic effect occurring at finite nucleation 
barriers. 
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I. INTRODUCTION 



Crystal growth takes place at the interface of the solid and liquid phases and is deter- 
mined by the structure of the crystal surface, the temperature and the supersaturation. 
Each growth unit is incorporated into the crystal when the energy barrier for absorption 
is overcome by the unit's kinetic energy. For rough surfaces, such as K- and S-faces, the 
density of privileged sites, i.e. kinks, is high and growth is fast. In this work, however, we 
consider initially flat faces such as the F-faces. For these surfaces, at low supersaturation 
values, kink sites are scarce and the growth process is slow. Growth units encounter the 
kinks directly from solution, or via a random walk due to 2D-diffusion on the surface. New 
steps are either generated through spiral dislocations or 2D-nucleation[l-6]. The tempera- 
ture is a determining factor. It is well known that there exists a critical value, Tr, where 
the step edge free energy vanishes and the crystal surface becomes rough at the growth unit 
scale. This is a thermodynamic phase transition, known as thermal roughening [7H1 2 1 . 

At constant temperature, the driving force is the chemical potential difference or su- 
persaturation and an increase in the driving force will translate into a decrease in critical 
nucleus size. Since for temperatures less than Tr th e step edge's free energy is greater than 
zero, according to the Gibbs-Thomson relation [I3I . Il^j ]. an activation barrier for the forma- 
tion of the critical two-dimensional nucleus exists and with it, a critical nucleus size larger 
than the crystal's individual building blocks. The size of this 2D-critical nucleus is inversely 
proportional to the temperature and the supersaturation |l[ . 
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16] that for temper- 



Recently, it has been reported for protein crystallization systems 
atures lower than Tr there exists a critical supersaturation for which the size of the critical 
two-dimensional nucleus is reduced to the order of one growth unit and the activation bar- 
rier for two-dimensional nucleation essentially vanishes. Due to a large step density and a 
very small two-dimensional critical nucleus at elevated supersaturations, the surface becomes 
rough and offers many favorable sites almost uniformly distributed across the surface. Con- 
sequently, arriving molecules can be incorporated quasi at any site. This transition from a 
slow, layer-by-layer growth regime to a fast continuous growth regime at high driving forces 
is referred to as kinetic roughening transition 
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22]. 



Kinetic roughening for the case of crystallization from solution has been observed for 
many small molecules, i.e. Si02, AI2O3, ZnO and ZnS 23|, |24j, n-paraffins 2M, [26[ and 
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many others. For the case of macromolecules such as in protein 



crysta 



lization, kinetic 



data of continuous growth has only been presented for lysozyme [15|, [27|, |28|. Very recently, 
quantitative data for both the classical layer-by-layer growth mechanism and the kinetic 
roughening regime for the crystallization of glucose isomerase from Streptomyces rubiginosus 
have been presented [it] ]. 

In this paper we show that kinetic roughening-like transitions in the presence of a fi- 
nite nucleation barrier can occur owing to kinetic mechanisms still rendering possible the 
transition from slow layer-by-layer growth regime (low supersaturations) to a fast continu- 
ous growth regime (high supersaturations). The general scheme used views the process of 
crystal growth at constant temperature driven by the difference in chemical potential of the 
solution around the crystal surface and the reservoir of the surrounding liquid phase. The 
local environment around the surface incorporates the barriers, kinks and dislocations while 
relevant kinetic processes such as absorption, surface diffusion and nucleation contribute to 
the growth. A mean field model and a Monte Carlo simulation are proposed and serve as 
an illustrative generic explanation of the observed data from protein crystallization 16]. 



Large reservoir 



Local environment 
around surface 



desorption 



Surface: absorption, surface diffusion, 
nucleation, growth 



In Section II a minimal mean field model is developed incorporating cluster formation 
and possessing the two main asymptotic regimes, exhibiting the transition from low to high 
supersaturation growth. In Section III we investigate an extension of a "solid to solid" 
(SOS) Monte Carlo simulation linked to, but more general than, the mean field model, and 
accounting for kinetic roughening in the presence of finite nucleation barrier and dislocations. 
The main conclusions are summarized in Section IV. 



II. MEAN FIELD MODEL 



Our mean-field model of crystal growth is based on a generalization of the Burton- 
Cabrera- Frank model for epitaxial growth [3, [29, [3^ and is also motivated by theoretical 
work hinting at a competition between surface aggregation and the formation of cluster 
by direct incorporation of material from solution|31]. Let Xq and X denote, respectively, 
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molecules in the bulk solution far from the growing crystal and in solution near the crystal 
surface; and Y, Cn the concentration of molecules adsorbed onto the surface but not incor- 
porated into the solid and of clusters constituting islands of solid material of size N. There 
are four basic physical processes that we consider. First is the diffusion of molecules be- 
tween the bulk solution and the solution near the growing crystal, X ^ X. Distinguishing 
between the bulk solution and that near the surface allows us to model the formation of a 
depletion layer as molecules are adsorped onto the surface and incorporated into the solid. 
The second process is that of adsorption from the solution onto the surface, X Y . We 
consider two pathways for the formation of an island of size N: adsorption and aggregation, 
X + {N — 1) Y Cn, and aggregation alone, NY <=^ Cn- Denoting the concentrations of 
X and X by x and Xq respectively, the surface concentration of adsorbates and iV-islands 
by y and c n , and assigning rate constants to each process gives the rate equations 

— = axo + k' y — (a + k ) x + k[cN — kixy N ^ (1) 

til/ 



^ = k x - k' y -(N-l) k lX y N - 1 - Nk 2 y N + ((N - 1) k[ + Nk' 2 ) c N 

d c N _ u „_.JV-1 , ,„ „.jV 

dt 



k\xy + k 2 y — k[cN — k' 2 CN 



The rate constants a and a' control the rate of exchange between the bulk and the 
solution near the crystal surface (in effect, this represents a zero- dimensional model for 
spatial diffusion). The constants ko and k' Q control the rates of adsorption and desorption. 
The remaining rate constants pertain to the rates of nucleation and dissolution of clusters. 
Note that the probability to bring together N molecules is taken to be proportional to the 
N-th power of the concentration: this is equivalent to assuming the rate goes like e~P NAW 
where AW is the work of formation per molecule of a cluster of size iV and using the low- 
density approximation that the free energy of a molecule in solution is proportional to the 
log of the concentration. 

Notice that in the perspective of Eqs.Q growth arises from the fact that the number 
of clusters Cn (N given) increases until the whole surface is covered. Thus, the total solid 
mass is M = Ncn and the growth velocity is dM/dt = dNcN/dt. (In this zero-dimensional 
model, there is no distinction between growth within a layer and growth perpindicular to 
the surface.) In actual fact, in addition to this mechanism one expects that the formation 
of bigger clusters from smaller ones through, e.g., the process Cn +Y—> Cn+i should 
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also contribute to growth. In this Section, we will limit ourselves to the model of Eqs.flT]), 
which will be used as a reference for sorting out the different growth regimes in a transparent 
manner and as a basis for comparison with the Monte Carlo simulations. This approximation 
is expected to be valid for high supersaturations where step growth is less important. Finally, 
we only consider the process in which one molecule is directly incorporated from the solution 
since, for the conditions considered here, we expect this to be dominant over processes 
involving two or more molecules coming from solution. 



A. A simplified model 

To illustrate the basic mechanism leading to a change in growth regimes, we first ignore 
depletion of the protein in solution near the interface so that a = a' = and x = xq and 
neglect all dissolution processes so that k' = k[ = k' 2 = 0. This limit should describe to 
a good approximation crystals growing in the kinetically-limited regime. What remains is 
then the simplified dynamics, 

^ = k x — (N — 1) k\xy N ~ l - Nk 2 y N (2) 

Since there is no detailed balance, there is no equilibrium state. However, a steady state 
where ^ = exists, implying a constant growth velocity, N^- = v. The concentration of 
adsorbed molecules is then determined by 

= k x - (N - 1) k lX y N - 1 - Nk 2 y N (3) 



and the velocity is 



v = Nk x xy N - 1 + Nk 2 y N (4) 



When the concentration of protein in solution is small compared to that adsorbed on the 
surface, x « y, one finds 

v = k x + ... 
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while the reverse circumstance, x » y, gives 



n , 

fc a; + ... 



iV 

There is therefore a crossover from small x behavior, where the aggregation process, Ny ^ 
cat, dominates with velocity v = kox to the large x behavior where the x + (N — 1) y cn 
process dominates with velocity v = j^zik x. Since this transition occurs in the normal 
growth regime - below the usual roughening transition - the growth rates are, in these 
limits, linear functions of the concentration, as expected. Small and medium sized clusters 
are those contributing the most in the switching, since for large N the growth velocities of 
the two regimes become practically indistinguishable. As we will discuss elsewhere, to get 
a more detailed sense of the interplay between these mechanisms, a simple model can be 
formulated for which numerically-exact results can be obtained using the theory of finite 



Markov processes 
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B. Full model 



In the simp 



investigations 
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ified model, all dissolution processes were neglected. Experimental 



34| have shown, however, that this is not valid for protein crystals growing 
at low driving forces. Here, we therefore incorporate desorption and dissolution processes in 
the full model. Including the desorption of surface molecules, k' ^ 0, does not change the 
picture dramatically: the small x behavior is no longer linear but there is still a crossover 
between one well-defined growth mode at small x to another, faster mode at large x as in 
the simple model. Including all of the evaporation processes, k[ 2 ^ 0, however, does change 
the picture significantly. First, detailed balance defines an equilibrium state at some finite 
value of Xq. For bulk concentrations above this value, growth occurs but it will eventually 
end when a steady (equilibrium) state is reached. In between, a quasi-steady growth regime 
occurs during which a growth velocity can be defined. Detailed numerical solutions of the 
model again show a crossover in the growth velocity as the bulk concentration increases. The 
requirement of detailed balance at equilibrium, i.e. that the forward rates and backward 
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rates of all processes are equal, gives 



axo = a'x* 



(7) 




kix y 



* *JV— 1 



hy 



*N 




where x* is the equilibrium value of x, etc. Given a value for the concentration in the bulk 
solution, these relations determine the equilibrium values of all of the other quantities as 
well as imposing the consistency relation 



Supersaturation at the surface is defined as log(x/x*). To illustrate the general behavior of 
this model, we have solved for the time-dependence of the various quantities for the case that 
a = a' = 0.04fc , h = 0.02fc , k[ = k' 2 = 0.004£;o, k' = 0.5k , k 2 = 0.01fc , N = 2 with initial 
conditions xq(0) = 1, and x(0) = y(0) = c(0) = 0. Here the original time t is normalized 
to k t. The value of k' reflects the idea that once adsorbed, a molecule spends on average 
a long time on the surface before being desorbed back to the bulk. The choice k 2 < ki 
amounts to assuming that the activation barrier for aggregation alone is larger than that for 
adsorption and aggregation. Aside from these assumptions and the physical requirements 
that the rates at which molecules join a cluster, k\ and are greater than the rates for 
the reverse processes and the detailed balance condition, EqJHJ these values are somewhat 
flexible. Given an initially flat surface, ctv(0) = 0, islands will form causing the mass of 
the surface to increase until, at long times, equilibrium is reached. At intermediate times, a 
more-or-less steady growth regime is reached, see Fig. [TJ corresponding to a plateau in the 
growth velocity, which defines a quasi-steady state growth velocity as shown in Fig. [5J 
Taking the maximum growth velocity as a measure of the typical velocity at intermediate 
times, one again finds a crossover between small x and large x regimes as illustrated in 



k' kik' 2 = kok[k2- 



(8) 



Figure |3j 
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FIG. 2: The dimensionless growth velocity, , as a function of time for the same parameters 
as used in Fig. [TJ The supersaturation in this case is log(x/x*) « —0.84. 



III. ONE-DIMENSIONAL KINETIC MONTE CARLO MODEL 



In order to study both the possibility of multiple growth regimes without vanishing nucle- 
ation barrier and the roughening transition, we have performed Kinetic Monte Carlo simu- 
lations of the growth of a simple, one-dimensional surface in the spirit of the Solid-On-Solid 
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FIG. 3: Maximal growth velocity as a function of supersaturation, log(x/x*), for the parameters 
given in the text. The broken lines extrapolate the values into asymptotic regions. 

(SOS) model|35[]. The system consists of a set of growth sites arranged as a one- dimensional 
lattice. There are N sites, each characterized by a height hi where i — 0, N — 1. Periodic 
boundaries are used so that formally, one has = ho. The idea is that molecules randomly 
land on the surface at some specified rate Vhit representing the diffusion of molecules onto 
the surface. Once on the surface, a molecule attempts to leave the surface at a rate u evap . 
Physically, u^u will be determined by the concentration of the molecules in solution and the 
rate at which they move while v evap is a measure of how fast a molecule makes an attempt 
to leave the surface. 

Nucleation is anomalous in one dimension due to the fact that the "surface area" of a 
cluster does not depend on its size. We therefore introduce a rule designed to capture the 
most important feature of surface nucleation which is the existence of a critical cluster size. 
Specifically, molecules are not allowed to leave the surface if they are part of a cluster that 
is of a specified critical size or larger. The critical size n c is a parameter. Molecule j is 
"part of a cluster" of size n if there are n contiguous sites, including site j, with heights the 
same as or greater than hj. Thus, small clusters can form spontaneously and dissociate but 
a cluster at or above the critical size is stable. 

A time step dt = 1/ {u hit + u evap ) is defined and the dynamics of the kinetic Monte Carlo 
algorithm are as follows: 
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1. At time t the system has some configuration of heights and total mass M (t) = ^ i=0 hi 
and roughness R(t) = X^o* \hi — Roughness could alternatively be charac- 
terized by the variance of the heights, the correlation length, ... 

2. The following is performed N times: 

(a) Choose a random site, j. 

(b) choose a random number u G [0, 1] and set hj — >■ hj + 1 if u < v hit dt. 

(c) if the height is not increased, (this happens with probability 1— u hit dt = u evap dt), 
then remove a molecule, hj — > h j — 1 if site j is not part of a supercritical cluster. 

3. Set t ->■ t + dt. 

4. Return to step 1 until the desired number of cycles is completed. 

As it stands, there are only two meaningful parameters: the ratio of the two rates, Vhit/v evap 
and the size of the critical cluster, n c . Since super-critical clusters are absolutely stable there 
is no equilibrium state: clusters form no matter what value the parameters are given, as in 
the simplified form of the mean-field model discussed above. The growth rate is calculated as 
(M (t + dt) — M(t)) /dt and is normalized to the number of sites to give the corresponding 
growth velocity. In the simulations, we also allow for the presence of "defects", localized 
regions in which the critical cluster size is smaller than elsewhere, to serve as sites for 
heterogeneous nucleation. 

Simulations have been performed under three circumstances: no defects, a "wall" defect 
and a "spiral" defect. The "wall" refers to a set of sites, say sites 0, n c for which the evap- 
oration rate is zero. These therefore grow very fast and serve as a source for heterogeneous 
nucleation. For the "spiral" defect, the critical nucleus for site is set to one, for site 1 it is 
set to 1, and so forth up to site n c . This is not meant to realistically model a spiral defect, 
but rather to test whether the results are sensitive to the shape and nature of the defect. 
The name indicates the similarity in shape of the resulting defect to a projection of a spiral 
defect onto two dimensions. In fact, we observe no qualitative difference in the results using 
the two different defects. The main effect is that at low supersaturation, Vhit/v e vap « 1, 
the rate of growth in the system with no defects is dominated by the time taken for nucle- 
ation to occur - by comparison, step growth happens relatively quickly. Thus, the measured 
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FIG. 4: Growth velocity 
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(panel a) and roughness (panel b) as functions of the hit rate for 



various choices of the critical size, n c in the presence of a spiral defect. 



rate of increase of the solid is basically a measure of the nucleation rate and is not directly 
relevant to the rate of step growth. However, the same qualitative effects are seen with and 
without defects. Intuitively, one expects that if v^t « v evav then nucleation will be rare 
and growth will be dominated by attachment to existing super-critical clusters, i.e. smooth 
step growth. When the hit rate increases, nucleation becomes more probable until at some 
point nucleation can occur easily and growth will be a combination of growth of existing 
clusters and nucleation of new clusters. Simulations were performed for a system consisting 
of 1500 sites and the velocity, denned as the rate of change of the total mass, was averaged 
over 1000 cycles. 

Figure H] shows the rate of growth of the surface as a function of the hit rate for various 
choices of the size of the critical cluster when a spiral defect is present. For a small critical 
cluster, n c = 2, nucleation is almost instantaneous for all but the lowest values of Vhujv evav : 
growth is dominated by nucleation and the growth velocity increases more or less linearly 
with the hit rate. For larger critical sizes, two different growth regimes can be distinguished 
in the Figure. At small hit rates, growth is dominated by heterogeneous nucleation at 
the defect with the subsequent growth of steps away from the defect. This growth by 
aggregation is analogous to the NY — > CV scenario in the mean-field model. However, 
at high hit rates, nucleation becomes important and the growth rate increases due to the 
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FIG. 5: Growth velocity as a function of hit rate for various choices of the critical size, n c with no 
defect. The initial growth rate is very small and is dominated by the nucleation rate. 



combination of processes with nucleation dominating at high hit rates. This is analogous 
to the (N — 1)Y + X — > Cn process in the mean- field model. The point of crossover 
between the regimes increases with increasing critical size and corresponds to a roughening 
transition as is illustrated in Fig. Hb- Note that for the smallest critical cluster size, the 
transition region is too small to see in this figure. For critical clusters larger than 2 and for 
small supersaturations, the roughness is virtually independent of the critical cluster size and 
increases slowly with the supersaturation while at higher supersaturation, the transition to 
the rapid growth regime is also signaled by a dramatic increase in roughness. 

Finally, Fig. [5] shows the growth rate when no defect is present. The main difference 
from the previous results is that at low supersaturations, the growth rate is very small since 
growth only occurs after the nucleation of a super-critical cluster. Once such a cluster forms, 
the growth of a new layer of crystal is quite rapid so that the overall growth rate is dominated 
by the nucleation rate. Nevertheless, a clear crossover from this slow growth regime, where 
most mass is added by step growth, to a fast regime, where nucleation plays an important 
role in adding mass, is evident. 
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IV. CONCLUSIONS 



Within the context of classical nucleation theory, kinetic roughening is understood as 
occurring when the supersaturation is so high that the size of the critical cluster becomes 
smaller than one molecule[l|. At low supersaturations, growth is dominated by the hetero- 
geneous nucleation of steps that then grow smoothly to cover the crystal surface while at 
high supersaturation, above the transition, homogeneous nucleation also occurs leading to 
rapid growth of a rough surface. In this paper, the possibility of a similar transition for 
finite-sized critical clusters has been investigated. Two models for the growth of crystal 
surfaces which show a transition from a slow to a fast growth regime were presented. In 
the mean field model, adsorption onto the surface with subsequent surface diffusion and 
aggregation competes with direct formation of clusters to give two different growth regimes 
for the formation of islands of new crystal. At low bulk concentrations, the adsorption and 
aggregation mechanism are dominant whereas at high bulk concentrations, the formation 
of critical clusters by the formation of subcritical aggregates with direct incorporation of 
material from solution dominates. The difference in growth rates is greatest when the crit- 
ical cluster is small (being a factor of 2 for a critical size of 2). These results elaborate on 
previous theoretical work that indicated the importance of direct incorporation of material 



from the bulk 
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A one- dimensional Kinetic Monte Carlo model based on the well-known SOS model 
showed both different growth regimes and roughening. The main difference from the usual 
SOS model is the definition of a critical cluster size to simulate nucleation as it actually oc- 
curs on a two-dimensional surface. This model showed a clear cross-over between a smooth, 
step-growth dominated regime at low supersaturation and a rough, nucleation-dominated 
regime at high supersaturation thus exhibiting the observed features of the kinetic rough- 
ening transition. The growth velocities both below and above the transition are roughly 
linear functions of the hit rate (Fig. H|), or in physical terms, the concentration, just as in 
the mean- field model (see, e.g., Eqs.flHl E])). 

The main conclusion from this work is that the kinetic roughening- like transitions do not 
necessarily require the vanishing of the critical cluster size. Instead, "kinetic" effects suffice 
to give qualitatively different growth regimes and roughening of the surface. In future work, 
it is intended to further elaborate the simple mean field model presented here to give a 
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semi-quantitative description of the results of simulation as well as of experiment. Finally, 
it is noteworthy that the ex per iments revealing the kinetic roughening mechanisms that 



motivated the present study [16|] were done using gel to quench convection which brings 
one close to microgravity conditions. It might be expected that a comparison under real 
microgravity conditions would reveal further important features. 
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